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Abstract 

We report the first observation of difltractive J/'iIj{-^ fi'^/i") production in pp collisions at y^=1.8 
TeV. Diffractive events are identified by their rapidity gap signature. In a sample of events with 
two muons of transverse momentum p^ > 2 GeV/c within the pseudorapidity region |?7| <1.0, the 
ratio of diffractive to total J/ip production rates is found to be Rj/^ = [1.45 ± 0.25]%. The ratio 
Rj/^{x) is presented as a function of s-Bjorken. By combining it with our previously measured 
corresponding ratio Rjj{x) for diffractive dijet production, we extract a value of 0.59 ± 0.15 for the 
gluon fraction of the diffractive structure function of the proton. 

PACS number(s): 13.85.Ni 

In the course of our studies of high energy pp interactions at the Fermilab Tevatron using the Collider Detector at 
Fermilab (CDF), we have observed a class of events incorporating both a hard (high transverse momentum) partonic 
scattering and the characteristic signature of single diffraction dissociation, namely a leading proton or antiproton and 
a forward rapidity gap, defined as the absence of particles in a forward pseudorapidity {rj) Q region. The rapidity gap 
in such "hard diffraction" processes is attributed to the exchange of a Pomeron g , which in QCD is a gluon/quark 
color-singlet construct with the quantum numbers of the vacuum. Experiments on hard diffraction can be used to 
address the question of whether the Pomeron has a unique, factorizable partonic structure. 

In four previous Letters, we reported results from diffractive W^-boson B, dijet M, and 6-quark pj production 
obtained using the rapidity gap signature, and dijet production in association with a leading antiproton l^]. These 
results include measurements of the gluon fraction |H,p| and of the xt,j (Bjorken x) dependence of the diffractive 
structure function of the antiproton |^ . Comparisons with diffractive deep inelastic scattering data obtained at the 
DESY ep collider HERA revealed a severe breakdown of QCD factorization, expressed mainly as a suppression of 
a factor of ^ 10 of the overall normalization of the diffractive structure function at the Tevatron. A breakdown of 
factorization was also observed between the diffractive structure functions measured from dijet production in single 
diffraction and in double Pomeron exchange at the Tevatron |7| . In this Letter, we report a measurement of diffractive 
J/^p production in pp collisions at y/s =1800 GeV, pp — *■ p(or p) + J/^ + X, and compare the J/tp diffractive fraction 
with the results of our previous hard diffraction measurements to further characterize the observed breakdown of 
factorization. 

The data used in this analysis were collected during 1994-95 and correspond to an integrated luminosity of 80 pb^^. 
The technique we use to extract the diffractive signal is identical to that used in our previous studies. In a data 
sample satisfying selection requirements for J/^ decaying into fi^ ii^ , we look for events with a rapidity gap in either 
of the two forward regions of the detector covering the pseudorapidity range 2.4 < \r]\ < 5.9. We define a rapidity 
gap as the absence of hits in the beam-beam counters (BBC), which cover the region 3.2< \ri\ <5.9, and the absence 
of calorimeter towers with energy above 1.5 GeV within 2.4 < |r/| < 4.2. The size of the calorimeter towers in this 
region is A?; x A^ = 0.1 x 5°. 

The detector components relevant to J/tp selection are the silicon vertex detector (SVX), the vertex time projec- 
tion chamber (VTX), the central tracking chamber (CTC), the central electromagnetic and hadronic calorimeters 
surrounding the CTC, and the central muon detectors. The J/^ acceptance is limited by the muon detectors, which 
cover the region \ri\ < 1.0. The SVX provides spatial measurements in the r-0 plane with a track impact parameter 
resolution of [13 -I- (40GeV/c)/pT] Mm- The VTX is used primarily to measure the longitudinal position z of an 
event's primary vertex, and the CTC provides momentum analysis for charged particles. The combined CTC/SVX 
transverse momentum resolution for charged particles is (Jp^/pT — 0.0009pT © 0.0066, where pt is in GeV/c. Details 
of the CDF detector components can be found in Ref. [^ . 

A three- level dimuon trigger system was used to select events with a pair of muon candidates M. At Level 1, 
the dimuon trigger selection required the presence of two radially aligned pairs of time-correlated hits in the muon 
chambers. Level 2 required that each pair of muon chamber hits match a track in the CTC found by the Central 
Fast Track (CFT) processor, which performed a partial reconstruction of all CTC tracks and determined the pr 
with a momentum resolution of (jp^lp\ ~ 0.035 (GeV/c)~^. At Level 3, performed in software, events were required 
to contain two oppositely charged muon candidates with an invariant mass within 300 MeV/c^ of the J/f/; mass of 
3096.9 MeV/c2 ]lO|. 



In addition to the J/f/) selection requirements used in previous CDF analyses |11|, the following two requirements 
were imposed on the data: first, since the BBC information is used to tag rapidity gaps, only data for which there 
was no BBC coincidence requirement in the trigger were considered; and second, since additional interactions in 



the same beam-beam crossing would most likely spoil a diffractive rapidity gap, only events with one reconstructed 
primary vertex were retained. In order to ensure that reconstructed muons were found in the kinematic region where 
the trigger is highly efficient, a minimum transverse momentum of 2 GeV/c was required for each muon candidate. 
For a precise vertex measurement, both muons were required to be reconstructed in the SVX detector. The dimuon 
invariant mass distribution for events passing the above requirements is shown in Figure la. The signal region, defined 
as the dimuon mass range of 3.05 < M^+^- < 3.15 GeV/c^, contains 18910 events. 

There are three sources of dimuons in the above J/ip candidate event sample: (a) J/ip^s directly produced in pp 
collisions, or resulting from decays of intermediate states which are sufficiently short-lived {e.g. from Xc decays) so that 
their measured decay vertex cannot be distinguished from the primary event vertex; (b) J/^'s from B-hadron decays, 
and (c) background from processes for which the dimuon invariant mass falls accidentally in the J/ip signal mass 
window; the latter includes dimuons from Drell-Yan production, double semileptonic 6-decays, meson decays-in-flight, 
and residual hadrons that penetrate the calorimeter and are misidentified as muons. The first two sources contribute 
to the J/ip invariant mass peak, while the last comprises the background under the peak. The fraction of background 
events in the signal region is evaluated by fitting the dimuon mass distribution with a sum of a Gaussian and a linear 
function. The fit yields a background fraction of (6.5 ± 0.1)% within the signal region. 

The fraction of J/t/i's from i?-hadron decays can be determined by fitting the proper decay length distribution, ct, 
using the appropriate function for each of the three dimuon components described above. However, because we do 
not fully reconstruct B decays, we use an approximation to ct described in [ [L1| and referred to as pseudo-cr. In the 
signal region, the fraction of background events is fixed at the value of 0.065, obtained from the dimuon mass fit, 
and the pseudo-cr distribution is fitted using for the background a parametrization derived from the sidebands and 
appropriate parametrizations for the prompt and B decay dimuon components P,ll| . The pseudo-cr distribution for 
the signal region is shown in Fig. lb with the fit result superimposed. The fraction of J/^p mesons from i3-hadron 
decays obtained from the fit is (16.8 ±0.4)%. The vertical line at 100 fan separates two regions: a "long-lived" region 
dominated by B decays, and a "short-lived" region mostly due to prompt J/^ mesons. The short-lived region contains 
15824 events, which are used in the analysis below. By numerically integrating the fitted B decay component in this 
region, the i?-hadron decay contamination is found to be 3.3%. 

As in our previous rapidity gap studies M-g, the diffractive signal is evaluated by considering the number of BBC 
hits, A^BBCi versus the number of the adjacent forward calorimeter towers with energy above 1.5 GeV, iVcAL- Figure 
2a shows the correlation between A^bbc smd Nqal- The multiplicity in this figure is for the side of the detector with 
the lower BBC hit multiplicity. The (0,0) bin, A^bbc = -^cal = 0, contains 92 events. The excess of events in this bin 
is attributed to diffractive production. The non-diffractive content of the (0,0) bin is evaluated from the diagonal of 
Fig. 2a with iVBBC = -^cal, shown in Fig. 2b. The non- J/i/; background in each bin of this plot, estimated by fitting 
the dimuon mass distribution to the sum of a Gaussian and a constant function, was subtracted from the number of 
J/ip candidates prior to plotting, yielding 87.4 ± 9.7 J/tp events in the (0,0) bin. An extrapolation to bin (0,0) of a 
linear fit to the data of bins (2,2) to (12,12) yields 19.9 ±3.9 non-diffractive events. The events in the (0,0) bin will be 
referred to as "diffractive" . Figures 2c and 2d show the J /if; transverse momentum and pseudorapidity distribution, 
respectively, for the diffractive and total event samples. In Fig. 2d the sign of the J/iJj pseudorapidity for events with 
a gap at positive rj is reversed, so that the gap appears always on the left. 

The number of diffractive events in the (0,0) bin must be corrected for the efficiency of requiring a single re- 
constructed primary vertex, ef^xi ^s well as for random BBC and forward calorimeter occupancy. The single-vertex 
requirement, which is used to reject events due to multiple interactions, also rejects single interaction events with extra 
vertices due to track reconstruction ambiguities. Removing the single-vertex requirement yields 15.4 more diffractive 
J/tp events, resulting in ef^^x — 0.85. For non-diffractive events, the efficiency of the single-vertex requirement, ei'Jj^xi 
was evaluated by comparing the ratio of single-vertex to all events with the ratio expected from the instantaneous 
luminosity. This comparison yielded eJIJtx = 0-56 ± 0.04. Finally, from a study of a sample of events with no re- 
constructed primary vertex collected in random beam-beam crossings, the combined BBC and forward calorimeter 
occupancy was measured to be 0.20 ± 0.06. 

After correcting the data for the efficiency of the single- vertex requirement and for the forward detector occupancy, 
we obtain a diffractive to total J/ip production ratio oi R^^F^ — (0.42±0.07)%. This ratio is based on diffractive events 
satisfying our rapidity gap definition. Therefore, it must be corrected for the rapidity gap acceptance, e^"^, defined 
as the ratio of events in bin (0,0) to the total number of diffractive events satisfying the same J/tp requirements and 
produced within a specified range of ^, where f is the fractional momentum loss of the leading (anti)proton. The 
gap acceptance for ^ < 0.1 was calculated using the POMPYT Monte Carlo generator ||l^ followed by a detector 
simulation. For a Pomeron structure function of the form /3/(/3) = 1//3 B, where /3 is the fraction of the momentum 
of the Pomeron carried by a parton, e^'^^ was found to be 0.29. Dividing i?j/^ by 6^"^ yields a diffractive to total 
production ratio of Rj/^ ~ (1.45 ± 0.25)%. 

The ratio Rj/^p is larger than the corresponding ratio for diffractive 6-quark production, i?^;, = (0.62 ±0.25)% g, by 



a factor of 2.34 ±0.35. As both J/i/; and 6-quark production are mainly sensitive to the gluon content of the Pomeron, 
we examine whether the difference in the two ratios could be attributed to the different average Xbj values of the two 
measurements. Given the x7°''^^ dependence of the diffractive structure function measured in dijet production [||, 

the double ratio Rt! = Rj/^/Rii, is expected to be equal to {xj- /x^)'^'^'^^ . Since in these measurements we 

consider only central J/t/j or 6-quark production, the ratio xj; /x^^ is approximately proportional to the ratio of 
the corresponding average pt value for each process, which is « 6 GeV/c for the J/-0 (see Fig. 2c) and « 36 GeV/c 
for the 5-quark (about three times the average pr of the 5-decay electron |^). The expected value for Rt! is then 
K, (6/36)^°""^ ~ 2.2, in agreement with the measured value of 2.34 ± 0.35. 

For a more direct study of the diffractive structure function, we restricted our analysis to events in which at least one 
jet was reconstructed. A jet is defined as a cluster of calorimeter towers within a cone size of Ai? = (A77^ + A0^)2 — 0.7 
with a seed tower of Et > 1 GeV. Since our diffractive J/^p events have a rapidity gap in the region 2.4 < \ri\ < 5.9, 
the core of the reconstructed jet for both diffractive and non-diffractive events is restricted to the region \ri\ < 1.7. 
The number of events passing this requirement is 8732. 

Figure 3 shows distributions for the J/ip + jet event sample: (a) is the diagonal (iVcAL, A^bbc) distribution, 
equivalent to that of Fig. lb, (b) the corrected S,p^p distribution for the (anti)proton on the side of the gap, evaluated 
using calorimeter and BBC information in a procedure described in Ref. H, (c) the J/ip transverse momentum, and 
(d) the azimuthal angle difference, 4> = \4>j/ip — 'Pjetl, between the J/^p and the highest Et jet. 

The Xbj of the parton in the (anti)proton participating in J/^ production is evaluated using the equation x^ = 

Py' (e^^-'/''' + e'^^^''*)/^/s, where the + (— ) sign stands for p{p). In leading order QCD calculations, the ratio 
of diffractive to non-diffractive production is equal to the ratio of the corresponding structure functions. For J/^p 
production, the ratio Rj/^{x) per unit ^ was evaluated for the events in the region 0.01 < ^ < 0.03 (see Fig. 3b) and is 
plotted in Fig. 4 along with the same ratio for dijet production, Rjj{x), obtained from Ref. [p|. The structure function 
relevant to dijet production is Fjj{x) ^ g{x) + ^q{x) Q, where g{x) and q{x) are the gluon and quark densities in the 
proton and „ is a color factor. For J/ip production, which is dominated by gg interactions, Rj/^{x) ~ g^{x)/g{x). 
The ratio of Rjj (x) to Rj/^ (x) is given by 

U f \ 1 + 4?fi£) 

^/^^ ' RjuiAx) 1+I2M ^ ^ 



^^/'a(-) 1+lfi 



where the superscript D is used to label the diffractive parton densities. Evaluating this ratio of ratios by integrating 
the Xbj distributions for Rjj and Rj/^ in the region 0.004 < a; < 0.01 (kinematic boundaries for full acceptance) yields 
[i?jj(x)/i?j/^(a;)]cxp = 1.17± 0.27(stat). Using this value in Eq. and the ratio of q{x)Jg{x) = 0.274 at x = 0.0063 
and Q = 6 GeV calculated from the proton GRV98LO parton distribution functions |1J] , the gluon fraction of the 
diffractive structure function of the (anti)proton is found to be /^ = 0.59±0.14(stat)±0.06(syst), where the systematic 
uncertainty includes in quadrature the uncertainties of all correction factors. This gluon fraction is consistent with 
the value 0.54 ± 0.15 obtained by combining the results of diffractive W, dijet, and &-quark production [g|. 
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FIG. 1. (a) Dimuon invariant mass, and (b) pseudo-cr for signal region with fit to the sum of contributions from prompt 
J /"i}) production, J/t/i mesons from B-hadron decays and non-J/tp background. 
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FIG. 2. J/?/> event sample distributions: (a) Beam-beam counter multiplicity, A'^bbc, for the BBC with the lower multiplicity, 
versus forward calorimeter tower multiplicity, A'^cal; (b) multiplicity distribution along the diagonal with A''bbc = NcAh in the 
plot of (a); (c) J/tp transverse momentum and (d) J/tp pseudorapidity. 
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FIG. 3. Distributions for the J/ip + jet event sample; (a) the diagonal of the A^bbc versus Ncal distribution with 
Nbbc ~ A^cal; (b) Pomeron beam momentum fraction, ^p^p (corrected), for events with the J/tp within \ri\ < 1.1 (the shaded 
area is the distribution for events satisfying the rapidity gap requirements); (c) J/^ transverse momentum; (d) azimuthal angle 
difference between the J/i/' and the leading jet. 
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FIG. 4. Ratios of difFractive to non-difFractive J/tp (circles) and dijet (triangles) rates per unit Cp(Cp) ^^ ^ function of 
a;-Bjorken of the struck parton of the p{p) adjacent to the rapidity gap. 



